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This article describes a model in which ferromagnetism necessarily accompanies a spin-density- 
wave lockin transition in the borocarbide structure provided the commensurate phase wave vector 
satisfies Q = {m/n)a* with m even and n odd. The results account for the magnetic properties of 
TbNi2B2C, and are also possibly relevant also for those of ErNi2B2C. 
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INTRODUCTION 



The material TbNi2B2C is one of a number of rare 
earth borocarbide materials that display a fascinating va- 
riety of magnetic and superconducting phases. Q This 
particular material has a phase transition to a spin- 
density- wave phase at Tn w 15 K followed by the a 
pearance of weak ferromagnetism at Twfm ~ 8 K 
A neutron diffraction study || of the magnetic struc- 
ture showed that the wave vector of the spin-density 
wave decreased with decreasing temperature in the spin- 
density wave phase until the transition temperature to 
the weak ferromagnetic phase was reached, at which 
point the spin-density-wave wave vector became approxi- 
mately constant at Q = .545a*. This behavior is charac- 
teristic of a continuous lockin transition to a commensu- 
rate phase occurring at Twfm, the commensurate wave 
vector being Q = .545a* = (6/ll)a*. This value for the 
locked in wave vector has been confirmed by recent high 
resolution magnetic X-ray scattering studies. Q] There is 
some hysteresis observed at Twfm, which is consistent 
with the idea that hysteresis can occur at a continuous 
commensurate-incommensurate transition (as well as at 
a first order commensurate- incommensurate transition). 

Similarly, in ErNi2B2C, there is a transition to a spin 
density wave phase, followed apparently by a transition 
to a weak ferromagnetic phase || as the temperature is 
lowered. (|, 0, ||, [| The study of the weak ferromagnetism 
is made more complicated in this case because the mag- 
netic transitions occur within a superconducting phase. 
Here also a lockin phase transition appears to accompany 
the transition to weak ferromagnetism. (The evidence 
for the lockin is that the spin-density-wave wave vector 
appears to become more or less independent of temper- 
ature below Twfm)- The lockin wave vector is close to 
Q = 0.548a* = (17/31)a* according to g and close to 
Q = 0.550a* = (ll/20)a* according to g^. 

A curious aspect of the behavior of both TbNi2B2C 
and ErNi2B2C is that the apparent incommensurate- 
commensurate transition and the transition to weak fer- 
romagnetism seem to occur at the same temperature. It 



is possible for two independent second-order phase tran- 
sitions to occur at nearly the same temperature if the 
interactions in the system happen to have just the right 
values so that this happens. However, this is very un- 
usual, and to have it happen in two different materi- 
als is even more unusual. This suggests that the lockin 
transition and the transition to weak ferromagnetism are 
not independent, but have a common origin. The pur- 
pose of this article is to describe a model in which weak 
ferromagnetism necessarily occurs simultaneously with a 
lockin transition in a spin-density-wave state, and that 
could therefore be relevant to the magnetic behavior of 
TbNi 2 B 2 C and ErNi 2 B 2 C. The model described here is 
related to that proposed in jnj . 

One of the results derived below is that weak ferromag- 
netism necessarily accompanies the lockin transition in 
the spin-density wave states of TbNi2B2C and ErNi2B2C 
only if the commensurate wave vector is of the form 
Q = (m/n)a*, with m an even integer and n an odd 
integer. An independent microscopic calculation of the 
nature of commensurate phases of ErNi2B2C, including 
the possible occurrence of ferromagnetism, has been car- 
ried out in JTT|for a number of different commensurate 
wave vectors. These results for specific wave vectors are 
consistent with the general rule concerning the occur- 
rence of ferromagnetism developed in this article. 

The work described below makes use of a Ginzburg- 
Landau type of analysis. This is a useful way of obtaining 
results of general validity which depend on the symmetry 
of the problem. Thus, the approach is complementary to 
the microscopic type of calculation described in JL1[. 



THE SPIN-DENSITY- WAVE PHASE 

A model specifically related to TbNi2B2C will now 
be described. Above the spin-density-wave transi- 
tion temperature the material TbNi2B2C has a body- 
centered tetragonal structure with space-group symme- 
try I4/mmm. The spin-density wave in TbNi2B2C is as- 
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sumed to be represented by the equation 

S(r)= [S„, l e mQi - r + S,^ e - m ^"-] ( 5(r-R-B). 

(!) 

Here R labels the unit cell, B = 0, ^(a + b + c) labels 
the magnetic ions in the crystallographic unit cell (at the 
corner and body center positions), n labels the harmonic, 
and i = a,b labels the symmetry-equivalent wave vectors 
Q a = (Q,0,0) and Q 6 = (0,Q,0), where Q « 0.55a*. 
Since the spin-density wave in TbNi2B2C is thought to 
be longitudinally polarized |l2), i.e. Sj is parallel to Qj 
(see however JD|), the nonzero components of the pri- 
mary order parameter Si^ (describing the first harmonic) 
are Si t a )X and Si t b, y , which will be denoted below simply 
by S x and S y , respectively. It is important to remem- 
ber, however, that S x and S y are complex numbers giv- 
ing the amplitude and phase of the spin-density waves 
with wave vectors Q a and Q^, respectively. (A simi- 
lar model relevant to ErNi 2 B 2 C would involve Si <a<y and 
Sifi,x since the spin-density wave in this material is trans- 
versely polarized. |(|) 

Given the above description of the spin-density wave, 
the Landau free energy describing a second-order phase 
transition to the spin-density-wave phase has the form 

F = A(\S x \ 2 + \S y \ 2 ) 

+B(\S x \ 2 + \S y \ 2 r + C\S x \ 2 \S y \ 2 , (2) 

where A cx (T — Tjv), with Tjv being the transition tem- 
perature. The constants B and C must satisfy B > and 
45 + C > for stability. If C < a tetragonal double- 
Q spin-density wave phase with |5a;| = \S V \ ^ occurs, 
while if C > an orthorhombic single-Q phase with ei- 
ther 1 5a; | ^ or I Sy I ^ 0, but not both, occurs. Below 
in TbNi 2 B 2 C, the structure becomes orthorhombic A Q 
so it is clearly the single-Q state that occurs in this ma- 
terial. A given material normally contains both single-Q 
domains (in spatially different regions of the crystal). In 
what follows we study explicitly only the properties of 
the single-Q domain characterized by a non-zero S x ; the 
properties of the domain characterized by S y follow im- 
mediately by rotating by 7r/2 about the c-axis. 

Spin-density- wave harmonics of order 2n + l with n an 
integer are induced by contributions to the free energy 
proportional to S 2n+1 , x (S*) 2n+1 + S% n+hx S% n+1 . Since 
this term is linear in 5 2n +i j2 :, the value of 5 2n +i,a- which 
minimizes the free energy must be non-zero. There are 
other terms in the free energy that are also linear in 
S2n+i,x and which also contribute to the induction of 
the (2n + l) th harmonic. However, since we are only in- 
terested here in demonstrating that the symmetry of the 
problem requires the existence of the (2n + l) th harmonic 
(and do not attempt to evaluate its magnitude) it is suf- 
ficient to consider only one example of the terms linear 
in S2n+i,x- Similarly, throughout this article, only one 



example of the type of term necessary for our purpose 
will be given. The wave vector associated with 5 2n +i !X 
is (2n + l)Q a - Similarly, charge-density-wave (or, equiv- 
alently, longitudinal lattice-displacement-wave) harmon- 
ics of order 2n and complex charge-density-wave ampli- 
tude pm are induced by terms in the free energy propor- 
tional to p2 n {S x S 2 x n ~ v Y +pl n S x Sl n ~ 1 . The conclusions 
of this paragraph are well known from studies of the spin- 
density- wave state of chromium. [|l5) 

LOCKIN TO WEAK FERROMAGNETISM 

Consider a single-Q spin-density-wave phase charac- 
terized by a nonzero S x , as described above. A lock-in 
transition is obtained by adding so-called lock-in terms 
proportional to 5£ and (5*) p (p is an integer) to this free 
energy. Since the free energy must be invariant with re- 
spect to a translation of the spin-density- wave by the dis- 
placement \ (a + b + c) , these terms are allowed only for 
of Q such that Q = 2(m/p)a* , where m (as well as p) is 
an integer. (The factor 2 in the expression Q = 2(m/p)a* 
appears because the unit cell is body centered.) Thus if 
Q is close to satisfying this commensurability condition 
for some m, it sometimes pays the system to to adjust 
its Q to exactly satisfy it, so as to be able gain energy 
from the presence of the lock-in terms in the free en- 
ergy. If p is odd, however, terms proportional to 5£ can 
not exist by themselves in the free energy since they are 
not invariant with respect to time reversal. The solution 
to this problem for odd p is to consider terms such as 
5 p 5oa;,|jl^| where Sq x is the ferromagnetic component of 
the spin density, which are invariant under time reversal, 
and since Q = 2(m/p)a* , are also invariant under body- 
centered Bravais lattice translations. Note that symme- 
try requires the ferromagnetic moment to be in the same 
direction as the spin-density-wave polarization vector [in 
this case both are directed along the x (or a) axis]. Thus, 
for p odd, the terms stabilizing a lockin transition are 

Fi oc kin = A Sq x + B \S x \ p cos(p(l))So x , (3) 

where A a > is expected, and S x = \S x \exp(i<f). These 
terms must be added to the original free energy given 
above. Because the second term in Fi oc ^ in is linear in 
Sq x , the free energy has its minimum at a non zero value 
of Sq x in the locked in phase, and weak ferromagnetism 
is thereby necessarily induced at an odd-p lockin phase 
transition. For a continuous lockin transition correspond- 
ing to even p, ferromagnetism will not automatically be 
induced, and thus would not be expected to occur simul- 
taneously with the lockin. It should be noted that the 
lockin terms of Eq. ^| exist in the commensurate phase 
free energy independently of whether the transition to 
the commensurate phase is first or second order. Thus, 
ferromagnetism necessarily occurs in any commensurate 
phase with Q = 2(m/p)a* where p is odd. 
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If i?o > is assumed, the minimum of Fi oc ki n occurs for 
p4> = 2m and Sq x < 0, or forp0 = (2r + l)ir and Sq x > 0, 
where r is an integer. Assuming p = 11, as is appropriate 
for TbNi 2 B 2 C (since Q = 2{m/p)a* = (6/ll)o* in the 
locked in phase), and taking r = 1,2, ...11 in pcf> = 2m, 
one finds eleven distinct commensurate domains, each 
with Sq x < 0. The new commensurate unit cell has a- 
axis length 11a, and the eleven disinct domains are ob- 
tained by starting with one of them, and then producing 
the others by translating the spin-density-wave structure 
by the original lattice constant a ten times. Translat- 
ing each of these different commensurate domains by the 
displacement \ (a + b + c) gives a corresponding domain 
with all of the spin directions (including the ferromag- 
netic component Sg x ) reversed in sign. An incommen- 
surate phase can be viewed as a sequence of commensu- 
rate domains separated by domain walls, which are often 
called discommensurations. fl6| In a continuous incom- 
mensurate to commensurate transition, the domain walls 
are swept out to the sample boundary as the temperature 
approaches the transition temperature. In a commensu- 
rate to incommensurate transition, domain walls must 
be nucleated in the commensurate phase to form the in- 
commensurate phase. Pinning of the domain walls, as 
well as the energy barrier required to nucleate a domain 
wall, contribute to hysteresis effects for both first order 
and second order lockin transitions. The lockin transi- 
tion will be first order if the interaction between domain 
walls is attractive, and second order if the interaction 
between domain walls is repulsive (e.g. as in p7[). 

Consider an incommensurate phase viewed as a peri- 
odic array of discommensurations, each of which sepa- 
rates two commensurate domains differing by a transla- 
tion of | (a + b + c) . Suppose also that the commensu- 
rate domains all correspond to a particular wave vector 
Q = (m/n)a* with m even and n odd so that each com- 
mensurate domain is ferromagnetic. However, from what 
was said above, neighboring commensurate domains have 
their ferromagnetic moments in opposite directions so 
that the incommensurate phase as a whole has no net 
ferromagnetic moment. Applying an external magnetic 
field in a direction parallel or antiparallel to the ferro- 
magnetic moments of the domains will then cause the 
domains with moments parallel to the external field to 
grow (by movement of the discommensurations) at the 
expense of those with moments antiparallel to the ex- 
terna field, thus giving a sort of field induced ferromag- 
netism to the incommensurate phase. This effect will be 
inhibited by discommensuration pinning. 

MAGNETIC EVEN-ORDER HARMONICS 

Now consider a magnetic state in which there is a fer- 
romagnetic moment Sq x and a spin-density-wave (which 
may be either commensurate or incommensurate for the 



purposes of this section) characterized by the complex 
amplitude S x . Spin-density- wave harmonics of (even) or- 
der 2n and complex amplitude S2n,x are induced in the 
combined ferromagnetic and spin-density-wave states by 
terms in the free energy given by 

Feven = Sq x [S2n,x (S x ) 2 ™ + S% n x S x n ] . (4) 

Since this contribution to the free energy is linear in 
S-2n,xi magnetic even-order harmonics are necessarily 
present in a state containing both a spin-density wave 
and ferromagnetism. Conversely, a state containing both 
a spin-density wave and one of its even harmonics, nec- 
essarily contains a ferromagnetic moment. 

It has been found that, in ErN^E^C, even-order mag- 
netic harmonics are present below JV-FM-(|, ^| This very 
unusual result is explained by the discussion of the pre- 
ceding paragraph. 

DISCUSSION 

The model described above appears to account well 
for the simultaneous continuous lockin and weakly ferro- 
magnetic phase transitions which occur at approximately 
8K in TbNi2B2C. In particular, the commensurate wave 
vector Q = (6/ll)a* is of the form Q = (m/n)a* with 
m even and n odd. Thus, the continuous lockin transi- 
tion is necessarily accompanied by the presence of weak 
ferromagnetism. The direction of the induced weak fer- 
romagnetic moment is required by symmetry to be in the 
same (or opposite) direction as the polarization vector of 
the spin-density wave, which is consistent with the ob- 
servation ofQ that the induced ferromagnetic moment is 
in the basal plane. If there are ferromagnetic domains 
present, the lowest free energy type of domain struc- 
ture is expected to have domains of different ferromag- 
netic spin orientations associated with commensurate re- 
gions shifted relative to one another by the displacement 
i (a + b + c) . Even harmonics have been observed very 
recently |L3| both above and below the transition temper- 
ature to weak ferromagnetism, and it would be of interest 
to determine the magnetic character of these satellites 
as a function of temperature. (In the model described 
above, even harmonics of magnetic character should oc- 
cur only in the weakly ferromagnetic phase). Also, ]lj| 
has presented evidence that the spin-density-wave phase 
of TbN^P^C may not be perfectly longitudinal, at least 
at some temperatures. A reorientation of the spin po- 
larization should not affect the general arguments of this 
article concerning the lockin to weak ferromagnetism. 

Some features of the weakly ferromagnetic phase ob- 
served in ErN^P^C are accounted for by the above work. 
For example, the even harmonics which occur in the fer- 
romagnetic phase |sj , and which have the unusual char- 
acteristic of being magnetic, are shown to necessarily ac- 
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company the presence of ferromagnetism combined with 
a spin-density wave. 

However, the observed magnitude Q of the locked in 
wave vector in ErNi 2 B 2 C does not obviously satisfy the 
requirement Q — (m/p)a*, with m even and p odd, de- 
rived above for continuous phase transition to a locked 
phase that is also ferromagnetic. The commensurate 
wave vector found in @ is Q = 0.548a* = (17/31)a*, 
whereas the commensurate wave vector found in Q |^] is 
Q = 0.550a* = (ll/20)a*. It should be noted, however, 
that the value Q = (28/51)a* = 0.549a*|ll|, and the 
value Q = (16/29)a* = 0.5517a* would both produce a 
locked in weakly ferromagnetic phase as described above. 
The question now arises as to what is the significance of 
variations of wave vector of the order of SQ « 0.001. The 
transition from an incommensurate to a commensurate 
phase is brought about by sweeping the domain walls be- 
tween different commensurate regions out of the crystal. 
If it turns out that the pinning of the domain walls does 
not allow this process to be completed, so that there re- 
main domain walls separated on average by a distance 
1000a, then this will cause a change in the observed Q 
of order 5Q ss 0.001a*. Clearly, the interpretation of the 
locked in phase of TbNi2B2C, where the commensurate 
wave vector involves smaller integers and the absolute 
precision of its measurement is less significant, is more 
straight forward than the interpretation of the locked in 
phase of ErNi 2 B 2 C. 

It is also of interest to compare the results of this arti- 
cle with those of the microscopic mean field calculation of 
|nj . Rcf. |ll| studied a number of commensurate spin- 
density-wave structures having Q close to Q — 0.55a* 
for a detailed model of ErNi 2 B 2 C. The commensurate 
structures of that article having Q = (6/ll)a* and 
Q = (4/7)a* both had ferromagnetic moments, whereas 
the structure with Q = (5/9)a* had no ferromagnetic 
moment. These results are consistent with the general 
rule stated above that weak ferromagnetism is necessar- 
ily present in a commensurate phase with Q of the form 
Q = (m/n)a* where m is even and n is odd. On the 
other hand, Ref. (ll[] finds a first order phase transition 
within the Q = (ll/20)a* commensurate phase from a 
non-ferromagnetic to a weakly ferromagnetic state. This 
does not contradict the general rule given above, since 
the rule does not require a commensurate phase of this 
wave vector to be ferromagnetic, nor does it prevent a 
particular choice of the model interactions from yielding 
ferromagnetism at this wave vector. In this case, how- 
ever, there is no reason for the ferromagnetism to occur 
simulataneously with the lockin (and |llj does not ar- 
gue that it does). Finally, it is interesting to note that 
certain aspects of the neutron scattering data of [|| were 
explained in pl| by invoking a commensurate phase with 
Q = (28/51)a*. Interestingly, the symmetry arguments 
of this article require that such a phase is weakly fer- 



romagnetic. (Ref. |T1J is silent on the question of the 
ferromagnetism of their Q — (28/51)a* phase.) 
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